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ABSTRACT 


The BOXPLX of a user specified cost function is evaluated 
as a control system design method. The emphasis is on compen- 
sator design for multiple-input-multiple-output plants, but the 
technique should be applicable to any control system with 
adjustable parameters. 

The engineer must define the plant dynamics, inputs and 
outputs, the desired output responses for the specified inputs, 
weighting factors for the performance measure, the tyve of com- 
pensation to be used, the free parameters, and initial values 
and constraints on the free parameters. 

The program gives a direct search of the optimum values, 
within the specified constraints, for the free parameters. 
Evaluation of the results is given in graphical plots of the 
output time responses including the desired output response as 
well as the compensated response. 

The Transfer Matrix Method is used to provide an analytical 
check on the accuracy of the method and the procedure is illus- 


trated with a two-input-two-output system example. 
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T EC INTRODUCTION 


As control systems become more and more complex, design 
engineers are faced with an increasingly difficult problem in 
the development of compensators.  Single-input-single-output 
linear systems can be compensated using Bode plots or root 
zes plots and trial and error methods to meet the system 
specifications [1]. When compensating multiple-input-multiple- 
Sueput linear systems, use »scanebe made of transfer matrix 
techniques described by Ogata [2] to achieve total decoupling 
and the desired output responses. Since decoupling during 
the transient period of real systems may require controls 
which exceed physical constraints imposed by the svstem, it is 
sometimes desirable to decouple during steady state conditions 
БЕС ow coupling during the transient period [3]. 

Non-linear, single-input, single-ouput svstems, wherein the 
non-linearity can be lumped into one element, can be analvzed 
and compensated effectively using the phase plane method (ror 


- 


second order systems) and describing function analysis 12,14) 


and become rather unwieldy, but do produce the desired results. 
In ail of the analysis and design approaches mentionec above, 
any one of several well established simulation programs can be 
used to evaluate the design. 

To reduce the tedium involved in the design of compensators, 
frequency domain and time domain computer optimization schemes 
have been develoved, which will set free varameters in the 


compensator to yield a system response which most nearly 








approximates a reference response specified by the user. Lima 
used this approach to address the problem of controller design 
for ships engaged in underway replenishment [5]. MacNamara 
applied this method to an autopilot design [6], and Vines 
developed a generalized program for compensator optimization 
in single-input-single-output systems [7]. 

The following chapters present a generalized discussion of 
transfer matrix analysis and ЗЕ СЕ Шын or multivariable 
systems, optimization techniques, and performance measures. 

A two-input-two-output system is then compensated using transzer 
matrix techniques, and using parameter optimization. The 


results are then compared. 





-- MUBDPVARIABLÓE SYSTEM COMPENSATION 


A. TRANSFER MATRIX-ANALYTICAL METHODS 
IN Inttroduetron 
A linear, multiple-input-multiple-output or multivariable 
system can be described by a transfer matrix, which is simply 
an extension of the transfer function concept. Consider the 


two-input-two-output open loop plant: 


RIIS) С1(6) 


R2(S] 9,5 


Figure 2-1 Multivariable Plant 


The Block Diagram of the plant can be represented in the 


general sense as shown in Figure 2-2. 


RIS) . Gen lS 6% ais 





R 2(S} 9,225] Б Cc 218) 


Figure 2-2 Generalized Multivariable Plant 





Algebraic equations for the outputs іп the independent variable 


S, can be written directly: 


Cl(s) eM ЫШЫ 23 Ерт2 (5) R2í(S) (I) 


C2 (s) 


Gp2l(s)Rl(s) + Gp22(s) R2(s) (2) 


If the inputs Rj(s) and the outputs C. (s) are each considered 
as vectors of one dimension (i, j = 1, 2) the equations can be 
written in matrix form: 

Cl(s) Gpll(s) Gpl2(s) Rl(s) 


C2(s) -Gp21(s) epee (sh) R2 (s) 


The number of inputs and outputs can be increased, and need not 


be equal so that in the general case: 


eus) Gpll(s) Gpl2(s) --- Gpij(s) Rl(s) 

e2(5) Gp21 (s) R2(s) 
1 = 1 (4) 
il ji 
à | jt 

С (а) Gpil(s) Gpij(s) Β. (5) 
ДИ + = J 

er: C (s) = GS) Rts) (5) 


where Gp{s) is the transfer matrix (ixj) of the plant. 


2. Feedback Compensation 
Transfer matrices can be used to describe the сспреп- 
sation of multivariable systems as well as the plant.  Consicer 


the two-input-two-ouput plant above with feedback compensation 


dd 








shown in Figure 2-3. 


HIS) 
ЖЕ) SO a 5 
Нә ©) 
Gp[S) 


Ἑπσαμς ο 
Multivariable Plant with Feedback Compensation 
Algebraic equations for R) (s) and R, (s) Can be written: 
(=) = ea e - Π12(5) C2(s) (6) 
БО 5) 5 02/5) и 21(5) С1(5) - 222(s) C2(s) (7) 
with C, R, and U two element, one dimension vectors these 


equations can be rewritten: 


Rl(s) Ul(s) H11(s) H12 (s) οι (5) 


Б2 (5) 92 (S) FALS) 522 (5) 2s) 





or in the general case: 


As) E US) (9) 
where H(s) is the transfer matrix (jxi) of the compensation. 
From the uncompensated plant: 
C(s) = G(s) R(s) (103 


EBESEt)tuting equation (9) for R(s) in Equation (10): 


C(s) = G. (s) l U(s) - H(s) C(s)] (11) 
Rearranging: 
[I * G, (s) H(s) ]|C(s) = б (5)0(5) ᾿ 


sce Gp(s) is ап іх) matrix, E(s) із а 1хі matrix the product, 


Gp(s)H(s), is a square, ixi matrix. Assuming that [I + Gp(s)H(s) ] 


m 


is non-singular the inverse can be taken. 


eea i plying both sides of equation™ (12) by [I+G (s)H(s)] Í; 
C (s) = [I + G (s)H(s)] IG (s)U(s) (13) 
C(s) = G(s)U(s) (14) 

The compensated system transfer matrix (ixj) is then: 

G(s) = [I + G,(s)H(s)] "G, (s) (15) 


The plant with feedback compensation can be depicted more 


clearly as shown in Figure 2-4. 








Figure 2-4 


Matrix/Vector Representation of Multivariable 

Plant with Feedback Compensation 
To determine the feedback compensation, H(s), required to 
achieve design specifications, the dynamics of the plant must 
be known and put into transfer matrix form, E ης The desired 
closed loop dynamics must then be reflected in the closed loop 
ransfer matrix, G(s). Then through matrix manipulation the 
compensation, H(s), can be computed in closed form provided 
that G(s) and G (s) are square and non-singular. 


Premultiplying both sides of equation (15) by [I+G(s) B(s) | 


and then postmultiplying by er 
1 + G,(s)H(s) = 6 (s)G(s)”” (16) 
Rearranging: 
SG, (s)H(s) - G,(s)G(s) ^ - І (17) 
Thenmultiplying both sides of equation (16) by G (s) Í 
EG) - Gs) T - G (7^ (18) 
3. Cascade Compensation 


Cascade compensation can be treated in a similar “fashion. 


Consider the same basic plant with cascade compensation as shown 





ШІ Еталге 2-5. 


ШЕ) ы 
C11 






216) 
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Figure 2-5 


Multivarıable Plant with Cascade Compensation 


The algebraic equations for R- (s) and R, (s) 


RI (S) 


Gell(s)U, (s) 


R2(s) Gc2l1(s)Ul(s) 


Spe Rl(s) Gcll(s) 
R2(s) Gc2l(s) 


in the general case: 


G(s) = Gp(s)Gc(s) 


C(s) = Gp(s)Gc(s)U(s 


+ бс12(5)02(5) 


522215 4425) 


Gcl2( 
Gee (s) үр 


) 


Gp[S) 


ares 


= 


ТЕ Sp(s) is a square non-singular matrix the compensation 


Gc(s) is given by: 


Gc(s) = Gp(s) G(s) 





The compensated plant can be depicted as shown in Базаке 2-6, 





Figure 2-6 


Matrix/Vector Representation of Multivariable 


Plant with Cascade Compensation 


4. Cascade and Feedback Compensation. 
Cascade and feedback compensation can also be combined 


SSNSHIOWD in Figure 2-7. 


EIS) (5) ФЕ) 


RIS 
<6, ч) A stg 


(5) 


Figure 2-7 


Matrix/Vector Representation of Multivariable 


Plant with Cascade and Feedback Compensation 


l6 








Us) = G_(s)E(s) (26) 


C(s) = G (s)U(s) = Gp(s)Gc(s) [R(s) = Н(5)С(з) | 6272 
C(s) = Gp(s)Ge(s)R(s) - Gp(s)Ge(s)H(s)C(s) (28) 
Gp(s)Gc(s) R(s) - [І J Gp(s)Ge(s)H(s)]C(s) (29) 


МЕЙЕР (5) is an ixj matrix Gc(s) will be a j x j matrix, and 


ПОЛА Бе jxi, in which case Gp(s)Ge(s)H(s) is a square 


matrix (ixi). If [I + Gp(s)Ge(s)His)] is non-singular, it can be 


inverted and: 


τ Gp(s)Ge(s)Hls)]"Sp(s)Ge(s)R (s) (ο; 


^ 


C(s) = G(s) R(s) 


The compensated system transfer matrix G(s) (ixj) is: 


G(s) = [I + Gp(s)GC(s)H(s) | : Gp (s) Ge(s) mum 
Equation (31) contains two unknown matrices, H(s) and Gc(s). 
One of these must be selected arbitrarily and then in certain 
special cases the other can be solved explicitly. For example: 


His) = |; i E (32) 


Imewnich case: 


G(s) = [1 + Gplsice(s)] " Gp(s)Cc(s) (33) 

[I + Gp(s)Ge(s)]G(s) = Gp(s)Gc(s) (34) 

Hc eos) ec Ss) ait = 52571 (ЗЭ) 

If Gp(s) is a square matrix G(s) and Ge(s) will be square. II 
E 


Gp(s) [G(s) 


- I], and G(s) are non-singular: 











Gp (s) ο S) s IJ (36) 


C зве (37) 


a 
а 
u 
Ι 
Q 
rÜ 
N 


If Gc(s) and (eis ^t - I]Gp(s) } are assumed to be non-singular 


equation (37) can be rearranged for convenience: 


Gp(s)Gc(s) - (cer ~ της (38) 
Gp(s) = [6 (s) + = lies Ge (s) 7 (39) 
сс(в)7% - [6(s)7? - 1] Gp(s) (40) 
Gc(s) = {{G(s)~+ - 1] ср(ѕ) 17: (41) 


This result can also be obtained directly from equation (37) 


using the matrix identity (ag) i = Dv а 


Бы m 


In each of the systems discussed above the compensation 


cannot be solved for in closed form, except in special cases 


as indicated. 








Б БАРАМЕТЕН ОРТІМІ2АТТОМ - СОМРОТЕН МЕТНОР 
1. Optimization Technique 

The оорусу Method of constrained optimization (BoxPLX) 
was utilized in this program. The algorithm suggested by 
Box (9) was modified by the programmer, R. R. Hilleary, Naval 
Postgraduate School, to find the constrained minimum, and 
includes an integer programming option, as suggested in (10). 
The Complex Method is a ШТ ОД efschesumeoensrnainedr, direct 
search, Simplex Method introduced іп (11). Direct search methods 
compare function and constraint values only in searching tor 
the minimum value of the function. These direct search methods 
have been widely used because of their robustness, reliability, 
and ease in programming and use. They have widesoread applic- 
ability. The one drawback is that they Es generally less 
efficient than gradiant-based techniques (12). A more detailed 
Se eıssion 0:2 the Complex Method is included in the program 
documentation. 

ЕП Піс огсагашп, the main program simulates the reference 
τας зала then calls FUNCTION BOXPLX which in turn calis FUNCTION 
Borer UNGTLON FE calls SUBROUTINE Plant which simulates the 
compensated system. FE then computes the performance measure 
and returns to BOXPLX. BOXPLA keeps track of the compensator 
parameter values, computes a new set of free parameters and 
EN c E again: This iteration continues until termination as 


described in the documentation. See Appendix A for more 


information on the program. 








2. Performance Measures 

In order for a numerical method tois parameter opti- 
mization to produce a result, a decision process must be com- 
pleted, the result of which is the set of parameters which is 
defined as "optimum." This term optimum is an enigima of sorts. 
The program addressed here utilizes a comparison of time domain 
response of a compensated system with a reference response 
specified by the user. The word optimum above is dependent on 
the designer's knowledge of what reference response is best 
med tO the output he is trying to control. Furthermore, 
with the exception of very simple linear systems, it is diffi- 
cult to match a reference response exactly, so the designer 
must decide 1f the comparison at steady state is more or less 
important than that during the transient period. Additionallv 
very high frequencies with periods less than one half the 
integration step size may be filtered out bv the digital process 
in the computer, and not be present in the simulated output. 
The designer is faced with some decisions which he must nake 
objectively before the "optimum" set of parameters can be 
computed. 

The decision process mentioned above is based upon а 

performance measure which is an indicator of the "goodness" 
of a given set of parameters being varied in the optimization 
process. The designer must select the performance measure; tne 
performance measure is then minimized by numerical methods, 
since the objective here is to have the system response match 


the reference response. 








The term selectivity is used to describe a quality of the 
performance measures. If three dimensional space is considered 
with the ο ποιος measure plotted as a function of two free 
variables in a compensator of the system, the resultant topology 
should include one or more minima for the performance measure. 
The minima with the lowest value are referred to as the global 
Minima, although the minimization performed is constrained. 
Selectivity is directly related c the steepness of the slope 
of the contour in the area immediately adjacent to the minima. 
This selectivity is a function of the performance measure used 
as well as the plant and compensator dynamics (2,8). 

Several performance measures have been suggested in numerous 
sources, some of which will be discussed here. In each case 
the performance measure will be represented by the variable J. 

Consider a single-input single-output system with output, 
c(t), and reference response, r(t). The error, e(t), for this 
case is defined as the difference between the reference res- 


ponse and the system response or: 
e(t) = r(t) - c(t) (42) 


M e eaa eerror the defining equation for inite 
time is: 
E 
J = e” (t)dt (43) 
O 


This performance measure is not very selective and tends to 


emphasize large errors and ignore small ones (2,8). 








το 5 Έτος mulciplied square-error. The defining 
equation for finite time is: 
E 
J = t e“(t)dt (44) 


Steady state error is penalized more heavily by this per- 
formance measure. This performance measure is said to be 
more selective than the integral square-error criterion in that 
the value of the integral tends to change more rapidly with 
changes in the system parameters (2,3). 

πο οι, οἱ tener or.. Tre defining equation for finite 
time is: 

J = Ге (ас (35) 
О 

This performance measure is slightly more selective than the 
integral square-error criterion (2,8). 


Кыыс о οἳ ειποοπο ειστε absolute error. The defining 


τη Гот Tinite time is: 


Ë 


J = Sec ga (46) 
О 


This performance measure is more selective than those mentioned 
previously and, as in the case with the integral of time multi- 
plied square-error criterion, this performance measure tends to 
emphasize errors late in the transient period and in steadv 


state and ceemphasize those which occur early in the transient 


period (2,8). 








The optimization of the performance measure can be addressed 
in two categories; the case where the reference response speci- 
fred is the ος Pines ion: or the reference response is some 
function, usually a second order response, which the designer 
believes to be a model response for his application. Іп the 
first case the peak overshoot, setting time, rise time, and 
natural frequency, are determined by the dynamics of the com- 
en ated system and by the Bee measure selected. How- 
ever, 1£ the designer desires a second order response with a 
certain minimal rise time, he can specify a second order 
reference response which meets this requirement, and use the 
integral absolute error or integral error-squared performance 
criteria to optimize the compensation parameters. 

Provisions for both of these approaches have been included 
in the program in that the user can reference a second order 
response, where he specifies the damping factor, natural frequency 
and gain, or he can utilize the table loop-up feature in which 
he can specify any reference response he desires. 

Other performance measures can be developed using classical 
measures of system performance such as maximum overshoot, time 
5-5 еггог tO reach Its first Zero, time to reach maximum 
overshoot, settling time, frequency of oscillation of the 
transient or steady state error. However, performance measures 
based on these characteristics would be complex and increase 
Semeucation time, and their desirability over those previously 


discussed is questionable. 


The performance measure used in the optimization program 








(see Appendix) is the sum of the weighted performance measures 


of each output of the multivariable system. 


The subscript i denotes the output for which Ji 


Jo P cul ated. 


The subscript k denotes one of four performance measures which 


can be selected оу the user to be applied to each output in 


CURN: 


Integral square error 
N 
Th ce ać 
me 1,1 
ісі 


Integral of time multiplied square error 


71,2 = 2 ari 12-28 ) 


Integral of absolute-error 
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Integral of time multiplied square-error 
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N represents the number of integration steps in the simulation, 


AT represents the integration step size, M is the number of 


outputs being optimized, and W is the weighting £actor applied 





to each output. 

It should be noted that the performance measures are an 
approximation and when the weighting factor is introduced 
the performance measure diverges significantly from the defining 
integral equation. Comparison of performance measures for a 
given plant with various types of compensation should be done 
Ehe same weighting factors. The same performance measure 
could be achieved for a system using two different compensation 
schemes, since the performance measures for each channel may 


paw significantly. 








Soe TWO ENPUT-TWO OUTPUT SYSTEM EXAMPLE 
ШӨ” Тгапстег Matrix Solution 

In order to demonstrate the program it was decided to 
analyze a simple system, synthesize the compensation required 
to produce a specified output, and then show that the optimi- 
zation program will arrive at the same results. It is emphasized 
that this program will not determine the type of compensation 
required to achieve the desired response. It will set free 
parameters in the compensators (specified by the user) to most 
nearly produce the desired response. 

The analysis and compensation was performed using the 
transfer matrix method described by OGATA (2). The uncompensated 
plant and the compensated plant were simulated to show that che 
compensation achieved the design specifications. Certain para- 
meters in the compensator were defined as free parameters. 
These free parameters were offset from their optimal values 
and constraints were introduced which restricted the range of 
values for these parameters during optimization. The plant 
and compensators with free parameters were simulated in the 
optimization program and optimization was accomplished. 

To illustrate, a simple two-input two-output linear 
system shown in Figure 3-1 was selected as the plant to be 


compensated. The design specifications chosen were: 


1. Channel one and channel two are to be decoupled during 
the transient period as well as in steady state. 
2. Channel one is to have a second order response to 


a step input with a natural frequency, N, m: 10, a 
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Figure 3-1 Multivariable Plant to be Compensated 
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τ гасгог, 5, of 0.4, and no steady state error. 

3. ἔΛΛ;; ΞΕ a second order response to a 
step input with a natural frequency, Wa’ of 4, a damping 
factor, C, of 0.6, and no steady state error. 

To achieve total decoupling the off-diagonal elements of 
the transfer matrix of the compensated system must be Zero. 


The desired compensated system transfer matrix, G(s), is: 
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10:0 0 
с2 + 85 + 100 
G(s) = 50) 
0 - 16 
L 57 + 4,85 + 15 


The input/output equations obtained from Figure 3-1 are: 
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КОШЕ = Шоп Of equation (53) into equation (54), and equation 


(54) into equation (53) produces equations (55) and (56). 
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Collecting like terms in equations (55) and (56) and rearranging: 
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Multiplying both sides of equations (57) and (58) Бу 
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The transfer Matrix for the uncompensated plant is: 


= (s+2) 2 
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The generalized compensation to be used is shown in Figure 
3-2. Equation (+1) will be used to calculate the required 
compensation; since in this case the feedback compensation 


Esonsser matrix, H(s), is the identity matrix. 
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(s) = {[S(s) - I]G pu mmn) 
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Since G's) is a 2x2 matrix in this example Ge (s) and G(s) will 


ΜΞ. ον inspection (Equation (52)), G(s) is non-singular 


and can be inverted. 








Тһе compensated system transfer matrix G(s) is given in 





11 (52). The inverse of G(s) is calculated by equation 
2). 
T 
G(s) р СОҒ (62) 
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Where cor? is the transpose of the cofactor matrix of G(s), 


and |G(s)| is the determinant of G(s). 
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The olant transfer matrix, G (s), πο ven in eauation (61) as: 
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PMúltiplying equation (61) by equation (66): 
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If this matrix is non-singular it can be inverted. 
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One of the possible realizations of this compensation is 
ӘНЕ ІП Figure 3-3. To obtain a confirmation of the design, 
the system was simulated with and without compensation using 
International Business Machines' Digital Simulation Language, 
ο ш Figures 3-4 through 3-9 show the results or that 
Semmeatlon, which indicate that the compensators used do in 
fact achieve the stated objectives of the design. In each olot, 

race number 1 is the desired response, Z, and trace number 2 
the system response, Y. In Figures 3-4B, 3-78, 3-8B, and 3-98 
the two traces are superimposed. In Figures 3-5A and 3-6A 
the desired responses are zero. The reader should note the 


scale factors present on some of these clots. 


r r amo r Optimization Solution 
The compensated system and the desired response curves 


were then simulated and solutions were obtained using the 

















Compensated Multıvariable Plant 


Figure 3-3 
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Figure 3-5A 


Y2 and Z2 vs Time with Rl = Unit Step, R2 = 0.0 
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Figure 3-583 


шапа ο. vs Time witi RI = Unit Step, 52 = 0.0 
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Figure 3-6A 


Yl and Zl vs Time with Rl - 0.0, R2 = Unit Step 
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Figure 3-7A 
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Figure 3-8A 


Yl and Zl vs Time with ВІ = Unit Step, R2 = ШЕЕ = Ер 


NK 


valence 21 


1.20 


1.00 
¡o ee 





0.80 
Mm A A ү НАНЫ 


0.60 


0. 40 


| AT V 
59.20 2.00 4.00 5.00 8.20 12.20 


Еее 5-5 


vs Time with Rl = Unit Step, R2 = Unit Step 
For Compensatec System 


38 





°ч 


(6222 
22 


0.40 0.60 


0.20 





ο 
3 
Өз 00 2.00 4.00 5.00 3.00 10.00 


είστε 
R2 = Unit Step 


v2 and Z2 vs Time with Rl - Unico r p, 


Y2 
gos 


0.80 


60 


0.40 





© 
ә 
sp. QQ 2.00 4.00 5.20 3.00 10.00 


Figure 3-9B 
v2 and 22 vs Time with Rl = Unit Step, 32 = Unit Step 
For Compensated System 


39 





program discussed in this section. Тһе time response results 
are shown in Figure 3-10. There is a slight difference between 
the desired ee curve and the output of the simulated 
system. 

Next, two parameter optimization was performed. The pro- 
portional gain constant of each compensator was selected as 
the adjustable parameters. Each of the four performance measures 
given by equations (47) through (51) were used in successive 
optimizations. All other aspects of the program were unchanged 
Па this optimization. Equal weighting was applied to each 
channel. The gains 100 and 16 were chosen because their 
variation should have similar impact on the two respective out- 
puts. The upper limits were 105 and 21 and the lower limits 
were 95 and ll respectively. The integration step size was 
.02 and the final time was 5 seconds. Optimization was begun 
with the two gains set at 105 and 11 respectively, and зі = R2 = 
Шы шер functions. The results are given in Table 3-1. 

The system was simulated using these values for the two 
free parameters. The simulation results are shown in Figures 
Aca rcouaa 3-14. This simulation indicates the time multiplied 
integral square-error performance measure produced tne best 


results. 
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Measure 
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Measure 





Integral 
Square- 
Error 


Time 
Multiplied 
Integral 
ралаге- 
Error 


Integral 
Absolute 
Error 
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Multiplied 
Integral 
Absolute 


Error 





1.2026097 x 1077 104.9999 20.99998 


IN MEE ET E 102.2554 15.27091 


EEES z 10 104.9862 1123227 


9.994698 x 10 ШЕ 2728 14.81198 


Table 3-1 
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Yl and Reference 1 vs Time After Gee nization With 
Integral-Square Error Performance Measure 


Cu 
O 








2 
22 


014 


012 


010 


008 


X10 


006 


004 


002 





“300 002 90% 906 008 gl 


ga güre 3=1.13 


Y2 and Reference 2 vs Time After Optimization Witn 
Integral Square-Error Performance Measure 
(Curve With Greatest Overshoot is 22) 





014 


012 


010 


008 


x107? 


006 


004 


002 


000 





200 002 sch 506 008 Т/с 


Eee L3 m br 


Yl and Reference 1 vs Time After Optimization With 
Time Multiplied Integral-Square Error Performance Measure 
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Y2 and Reference 2 vs Time After Optimization With 
Integral Absolute Error Performance Measure 
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CONCLUSIONS 


From the results of the two-input two-output example 
system, the following conclusions can be stated: 

1. The technique of cost function minimization in the 
time domain can be used effectively for compensator parameter 
optimization in multiple input multiple output control systems. 

2. Cost function minimization is not a substitute for 
control system design. The user must be able to select the 
proper type of compensator which is capable of achieving the 
desired results, before initiating the optimization or the 
free parameters. 

3. The four performance measures evaluated produced 
slightly different solutions to the same problem, with the 
time multiplied integral square-error yielding the kest 


results in the example given. 
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